#include <stddef.h>

void *solve(size_t n, 
		size_t m,
		size_t k,
		const float *a,
	        size_t lda,
		const float *b,
		size_t ldb,
		const float *c,
		size_t ldc);


#include <stdio.h>
#include <time.h>

#define N 128
#define M 128
#define K 128
#define aS 4
#define bS 4
#define cS 4

int main()
{
	int n = 33, m = 23, k = 29;
	float a[M][K], b[K][N], c[M][N];
	for(int i = 0; i < m; i ++)
	{
		for(int j = 0; j < k; j ++)
		{
			a[i][j] = j * i + 1;
		}
	}

	for(int i = 0; i < k; i ++)
	{
		for(int j = 0; j < n; j ++)
		{
			b[i][j] = j * i + 9;
		}
	}


	for(int i = 0; i < m; i ++)
	{
		for(int j = 0; j < n; j ++)
		{
			c[i][j] = 0;
		}
	}
	

	solve(n, m, k, a, aS, b, bS, c, cS);

	for(int i = 0; i < m; i ++)
	{
		for(int j = 0; j < n; j ++)
		{
	
			printf("%5f ", c[i][j]);
		}
		puts("");
	}	
	return 0;
}
